Angular momentum transfer torques in spin valves with perpendicular magnetization 
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Spin valves incorporating perpendicularly magnetized materials are promising structures for mem- 
ory elements and high-frequency generators. We report the angular dependence of the spin-transfer 
torque in spin valves with perpendicular equilibrium magnetization computed by first-principles 
circuit theory and compare results with experiments by W.H. Rippard c.s. [Phys. Rev. B 81, 
014426 (2010)] on the CoFe|Cu|CoNi system. Furthermore, we predict a non-monotonous ("wavy") 
spin-transfer torque when the Cu spacer is replaced by a Ru layer. 

PACS numbers: 72.25.Ba, 85.75.-d, 72.10.Bg 



I. INTRODUCTION 

A current can be used to read out the information 
in magnetic memory devices by the giant magnetoresis- 
tance. Magnetic random access memory (MRAM) tech- 
nology has become scalable by writing information us- 
ing the current-induced spin-transfer torques (STT)ii~— 
The critical electric current density j c necessary to switch 
a magnetic layer in a spin-valve structure is an all- 
important figure of merit in this case. The introduc- 
tion of materials with perpendicular magnetocrystallinc 
anisotropics that forces the equilibrium magnetization 
out of the plane r- has helped to reduce j c ^- — 

Co|Ni multilayers are an interesting system with per- 
pendicular anisotropy ) 6 i 9 i 10 with a higher polarization 
and less spin-flip scattering than, for example, CoPt 
alloyAi Rippard et alJ^ studied current-induced high- 
frequency generation in structures with a perpendicularly 
polarized (Co|Ni)„ multilayer serving as the switchable 
magnet and an in-plane magnetized Co layer as polar- 
izer. The output power of such a device depends sensi- 
tively on the asymmetry of the angular dependent STT 
when the magnetization of the free layer is reversed] 1 ! 2 ! 13 
By generating an rf output by a dc current in a spin 
valve in which the free layer is magnetized normal to the 
polarizing layer, Rippard c.s. parameterized the skew- 
ness of the torque as a function of magnetization angled 
Koyama et alJ^- measured high speed current-induced 
domain wall velocities (40 m / s) in magnetic perpendic- 
ular Co|Ni multilayers with current-in-plane configura- 
tion. Another interesting materials system is Co|Ru ) 15 i 16 
which also displays perpendicular magnetic anisotropyJ^ 

Scmiclassical theorie s 18 ' 19 that combine a quantum 
treatment of the interface scattering and diffusion treat- 
ment of bulk scattering in general explain experiments 
on magnetic metallic multilayers welli^ Here we report 
calculations of the STT of spin valves containing per- 
pendicularly oriented ferromagnetic materials based on 
magnetoclcctronics circuit theory using interface trans- 
port parameters computed by first principles. The spin- 



orbit coupling is the origin of the magnetic crystalline 
anisotropy and perpendicular magnetization. However, 
the experimental spin-dependent interface resistances for 
not too heavy elements can be reproduced by parameter- 
free calculations without taking into account the spin or- 
bit interaction^ which will therefore be disregarded in 
the following. 

Here, we study the angular dependent STT in CoiNLj 
(the subscripts refer to the number of atomic lay- 
ers) based spin valves by circuit theory in combi- 
nation with first-principles calculations. Firstly, we 
present results for Co|Cu|(CoiNi x )j,Coi|Cu(lll) stuc- 
tures, where the subscripts 1 and x indicate again the 
number of atomic layers, while y is the number of 
stacks and compare them with experiments^ 2 - Next, 
we report large and 'wavy' angular-dependent STT for 
Co|Ru|(CoiNi 2 )a;Coi|Ru(lll) spin valves which might 
therefore be very efficient high-frequency generators. 

In Sec. HH we introduce our method to calculate the 
STT in spin valves in terms of the spin mixing conduc- 
tances of the interfaces computed from first principle, 
including corrections for the magnetically active bulk 
material and the diffusive environment. In Sec. IIIII we 
present results for the spin mixing conductances for two 
the two types of spin vales with perpendicular magnetic 
anisotropy and compute the angular dependence of STT 
by magnetoelectronic circuit theory. We summarize our 
results in Sec. IIV1 



II. SPIN MIXING CONDUCTANCE IN A 
DIFFUSIVE ENVIRONMENT 

The STT due to a current bias / in 
ferromagnet|normal-metal|ferromagnet (F|N|F) spin 
valves in which the magnetizations are at an angle 
can be computed analytically by circuit theory^! and, 



assuming structural symmetry, be parameterized a: 
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where the asymmetry parameter can be expressed in 
terms of the the parameters of the N|F interface as 

A = (771/^(1 -P 2 ) Re 77, where f, = 2G n / (G t + G^j 

is the normalized effective spin-mixing conductance and 
P = — Gij I {Gf\ + Gij is the conductance polar- 
ization. Here, G^,G± and Gfi are the spin-dependent 
and spin-mixing conductances, respectively, where the 
tilde indicates that they have been "Schep corrected" 
for a diffusive environment and include the effects of the 
magnetically active contact regions close to the interface. 
In deriving Eq. ([T]) spin flip in the normal layer has 
been disregarded. When the spin- flip diffusion length in 
the magnetic layers is much longer than the bulk layer 
thickness:— 
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where a is the spin index, dp(jv) the thickness of ferro- 
magnet F or normal metal N layer, p the bulk resistivity 
(for a single spin), and Ap the pillar cross section. The 
G sh, s are Sharvin conductances, G-j- = (e 2 //i) trtlt-j-, 

G i = (e 2 /h) trtjtj. and G n = (e 2 /h) tr (i-rjr^), 

where tm) ( r t(4-)) are the matrices of the transmission 
(reflection) coefficients of the phase coherent region of 
the N|F contact as seen from the normal metal and at 
the Fermi energy. I is an M x M unit matrix, where M 
is the number of conducting channels in N. The third 
term on the right-hand side of the last two equations 
are the Schep correction, while the second terms correct 
for the magnetically active bulk regions. When the ferro- 
magnetic layer is much thicker than the spin-flip diffusion 
length lf d , the latter should replace dp in Eq. (J2|). With 
the spin-orbit interaction we also ignore intrinsic spin-flip 
scattering at the interfaces. The ferromagnetic layers are 
assumed sufficiently thick such that mixing transmission 
contribution may be disregarded^ Note that Eq. ([T]) 
only holds for structurally symmetric spin valves. In the 
following we use the general expression in which the left 
and right interface parameters differ, as shown in FigfTJ 
but do not list the expressions explicitly here (see Refs. 

In our calculations the atomic potentials were deter- 
mined in the framework of the tight-binding (TB) lin- 
ear muffin-tin- orbital (MTO) method^ based on den- 
sity functional theory in the local density approximation 
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FIG. 1: Scheme of asymmetric F1|N|F2|N spin valves with 
perpendicular magnetization F2 used in the calculations 



and an exchange-correlation potential parameterized by 
von Barth and Hedin.— The self-consistent crystal po- 
tentials were used as input to a TB-MTO wave-function- 
matching calculation, from which we obtained the trans- 
mission and reflection at the interfaces. The calculations 
are carried out with a k| | mesh density equivalent to more 
than 3600 k|| mesh points in the two dimensional Bril- 
louin zone (BZ) corresponding to the interface unit cell. 
The technical detail can be found in Ref. HE Table U 
compiles our results for various interface conductances 
including the bulk corrections due to magnetically active 
regions. 



III. COiNI 2 |CU AND CO1NI2 RU 
MULTILAYERS 



We first focus on the Co|Ni multilayers, which we 
treat as phase coherent regions, i.e. we compute the 
scattering matrix of the entire multilayers, which is 
then treated in the circuit theory of conventional spin 
valves just like a single interface. We present the spin- 
dependent and mixing conductances of Cu|XjCu with 
[X„= (Co 1 Ni 2 ) n Coi]. Here the Cu leads on both sides 
are semi-infinite. X n denotes n repetitions of the C01M2 
multilayer unit. As in the experiments,— a Co atomic 
layer is added for better contact with the Cu reservoirs. 
Since samples have been grown by sputtering, we take 
interface disorder into account, which is in general well 
modeled by a two monolayer 50%-50% interfacial al- 
loy (CoiNi 2 )„ -> ([Coo. 5 Nio. 5 ]Ni[Coo. 5 Nio. 5 ])„^ Spin-flip 
scattering at the Co|Ni interface will suppress any ben- 
efits of an even larger number of Co|Ni interfaces.— We 
therefore present here only calculations with n < 5. The 
computed dimcnsionless mixing conductance f] is also 
listed in the table. 

In the fee crystal structure Co and Ni have nearly iden- 
tical band structures for the majority spin, which results 
in very transparent Co|Ni interfaces. The majority spin 
conductance therefore stays nearly constant with increas- 
ing n. For minority spin electrons the scattering at the 



FIG. 2: Comparison of computed and experimental (Ref. Il2i) 
angular dependent STT in Fl|Cu|(CoiNi 2 )5Coi|Cu (Fl=Co, 
CogoFeio) spin valves with two monolayer 50%-50% inter- 
mixed interfaces. For Co as fixed lead, we vary the thick- 
ness dco from 5nm to 20 nm. When using CogoFeio as fixed 
layer, we use rfcoonFem = 2.5 nm, and resistivity pco ao Fe 10 = 
154f2nm (Ref. [291) and Co|Cu interface parameters. The 
dark area indicated the experiment results parameterized by 
Slonczewski's formula with A = 1.3 and A = 1.7. Calcu- 
lations are carried out by circuit theory for an asymmetric 
spin valve with first-principles interface parameters using the 
Schep correction including the contribution form the magnet- 
ically active region of the bulk ferromagnet as described in 
the text. 
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Co|Ni interface is much stronger. Consequently the mi- 
nority spin conductance decreases rapidly with increasing 
number of Co|Ni interfaces. 

Fig. [5] shows the angular dependent STT exerted on 
the right-hand side of Fl|Cu|(CoiNi 2 ) 5 Coi|Cu (Fl = Co, 
C090FC10) spin valves with intermixed interfaces calcu- 
lated by our first-principles circuit theory and compared 
with the experimental result^ For pure Co as fixed lead, 
we vary dc from 5 — 20 nm and find that the angu- 
lar dependent STT falls into the experimental ranged 
estimated by Slonczewski's formula for symmetric spin 
valves with A = 1.3 and 1.7. Experimentally, CogoFeio 
is used as fixed layer. Its spin-flip diffusion length is 
shorter than that of Co, but its resistivity is also higher, 
so there is not much difference when compared with a 
Co polarizer. We assume that the interface is not af- 
fected. We plot the results of CoFe in Fig. [2] with 
dco 90 Fc 10 = 2.5 nm (Ref. HU), and p C o 90 Fc 10 = 154 flam 
and is very similar to pure Co with dc = 5 nm. The 
results for Co|Cu|(CoiNi a; )j,Coi|Cu from x = 2 — 3 and 
y = 2 — 5 arc shown in Fig. [3l We observe large differ- 
ence between epitaxial and disordered samples, but only 
weak dependences on x and y. The results for epitaxial 
(disordered) samples fall into the range of Slonczewski's 
A = 1.05-1.15 (1.4-1.5). 

The experimental results were parameterized by Eq. 




60 90 120 

Angle (degree) 

FIG. 3: Angular dependent STT in Co|Cu|(CoiNi 2 ) I Coi |Cu 
and Co|Cu|(CoiNi3) x Coi |Cu spin valves with dco — 5 nm. 



(TTJ) for a structurally symmetric spin valve, whereas 
our results are based on the theory for asymmetric 
structures i^i We suggest that in future experiments Slon- 
czewski's formula should be replaced by a more accurate 
parameterization. 

Another interesting material with perpendicular mag- 
netic anisotropy is Co|Ru. Experimentally, both 
hcp(OOOl)^ and fcc(lll)^ structures have been re- 
ported. Despite the large lattice mismatch between Co 
and Ru, hep Co|Ru could be grown epitaxially and the 
magnetic anisotropy depends on the thickness of the Co 
layer^ However, the metastable structure relaxes to 
a more stable one after annealing^ Co|Ru|Co with a 
metastable fcc(lll) structure has also been reported.— 
Here we present systematic calculations of the transport 
properties of Co|Ru pillars with different structure and 
lattice constants as listed in Table HT1 

For epitaxial samples, we show results for an fcc(lll) 
texture with lattice parameters for Ru, Co, and its av- 
erage. The lattice parameter along the growth direction 
is varied to keep the atomic volume constant. Both spin 
polarization and specific resistance are close to the ex- 
perimental values,— but considering the large lattice dis- 
tortion (7.3 ~ 14%) this may be accidental. 

For the epitaxial hcp(OOOl) texture our calculations 
yield very high spin polarizations P = —39 ~ —55% 
for both clean and dirty interfaces when Co adopts the 
Ru structure and lattice constants as reported^ and 
small specific resistances AR = A/G = 0.69 — 0.78 x 
10 _15 Om 2 . Here and below P and G have been Schep 
corrected with magnetically active layer thickness dc = 
5 nm. Note that the structure is metastable and under 
annealing Co is expected to return to its normal lattice 
parameter. 

To simulate sputtering conditions, a 14 x 14 Co is 
matched to a 13 x 13 Ru lateral super-cell for both 
fcc(lll) and hcp(0001), leading to a spin polarization 
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FIG. 4: Angular-dependent torkance T/I on the right- 
side ferromagnet in disordered Ru based spin valves with 
X n = (CoiNi2)nCoi. The Co and Ru buffer layers are assumed 
much thicker than the spin diffusion lengths and i^T, so 
that the Schep correction includes the bulk scattering of the 
latter length scales. We disregard the bulk scattering in the 
Ru spacer, which should be allowed for the small thickness of 
8 monolayer (2.21 nm) considered. 

of P = —15% and specific resistance of AR = 0.75 x 
1CT 15 fim 2 for a clean fcc(lll) texture, P = -28% and 
AR = 0.93 x 10~ 15 fhn 2 for a clean hcp(0001) texture. 
A 50%-50% interface alloy has little effect on fcc(lll) 
texture, but leads to a reduced P = —19% for the 
hcp(OOOl) texture. The measured spin polarization for 
the Co|Ru interface is P = —20% with specific resistance 
AR = 0.5 x 10" 15 Qm 2 (Ref. ? ). 

In Table U we observe that in contrast to the Co|Cu 
interface, Co|Ru has a negative spin polarization for both 



fcc(lll) and hcp(OOOl) orientations. Interesting is the 
relatively large dimensionless mixing conductance rj. The 
predicted very large mixing conductance implies a large 
skewness of the angular dependent STT, which makes 
this material promising for applications in high-frequency 
generators. 

Fig. 2] gives the angular dependent T/I in 
Co|Ru|FM|Ru(lll) spin valves. Here, disorder is mod- 
eled again by two monolayers of a 50% — 50% inter- 
face alloy (to rather small effect) and Schep and mag- 
netic bulk corrections have been implemented. When 
fitted by Slonczewski's formula, the STT on the soft 
C01N12 multilayer in the strongly asymmetric spin valve 
Co|Ru| (C01N12) Coi|Ru(lll) shows a large variation in 
skewness in terms of the parameter A = 0.5 — 2.0. The 
maximum of the angular dependent spin torque is shifted 
gradually from low angle to high angle when the thick- 
ness of CoiNi2 increases from 2 to 6 periods. When Co 
serves as the free layer, the (modulus of the) angular 
dependent torkance shows two peaks and a compensa- 
tion point when the thickness of CoiNi2 increases from 
3 to 5 periods. This shape can be understood in terms 
of the spin accumulation in the normal metal spacer in 
the parallel configuration^ which is accompanied by a 
non-monotonic angular magnetoresistance. This behav- 
ior has been observed in Py|Cu|Co and dubbed 'wavy 
torques'. 41 ' 42 



IV. SUMMARY 

We studied the angular dependent STT for materi- 
als with magnetization normal to the interfaces by cir- 
cuit theory in combination with first-principles calcula- 
tions. An interesting angular dependent STT is found 
in Co|Ru|(CoiNi2)a;Coi|Ru(lll) spin valve. Moreover, 
a 'wavy' angular-dependent STT acts on the Co layer 
in Co|Ru|(CoiNi 2 )xCoi|Ru(lll) structures. When the 
CoNi is the free layer, we expect very efficient high- 
frequency generation. 
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TABLE I: Parameters for clean (disordered) interfaces (in units of 10 15 fl^ 1 m~ 2 ). X n = (Co 1 Ni 2 ) Coi;Y n = (Co 1 Ni 3 ) n ; fee 
Cu and Ru have Sharvin conductances of G^ u = 0.55 x 10 15 fi -1 m~ 2 and G^ (/cc) = 0.98 x 10 15 n _1 m~ 2 , respectively. 
Cco(/cc) t = an( i ^Co(fcc) i = 1-09 x ^0 15 m~ 2 for majority and minority spin in fee Co, respectively. For hex Ru's 
Gtu(he X ) = 0.80 x 10 15 fT 1 m~ 2 and Co with hex Ru structure G 3 c h o{hexht = 0.40 and Gco(/«»U = 0.80 x 10 15 O -1 m~ 2 for 
majority and minority spins, respectively. A magnetically active bulk region correction is implemented for the normalized spin 
polarization P and relative mixing conductance fj. For Cu|Co interface, we use dco = 5 nm (Ref. [13 ). and bulk resistivity 
PCo = 60 f2nm with spin asymmetry ft = 0.46, which results in pQ = 81 Slnm and Pq d = 219 finm (~Ref l3lh . We use a spin 
diffusion length = 60 nm. 
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TABLE II: Comparison of the calculated spin polarization 
P (we use dco = 5nm) and specific interface resistances 
AR — A/G of clean (disordered) Co|Ru for different lattice 
parameters with experiment. 
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